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We discuss a technique for solving the Landau-Zener (LZ) problem of finding the probability of 
excitation in a two-level system. The idea of time reversal for the Schrodinger equation is employed 
to obtain the state reached at the final time and hence the excitation probability. Using this 
method, which can reproduce the well-known expression for the LZ transition probability, we solve 
a variant of the LZ problem which involves waiting at the minimum gap for a time t m ; we find an 
exact expression for the excitation probability as a function of t w . We provide numerical results to 
support our analytical expressions. We then discuss the problem of waiting at the quantum critical 
point of a many-body system and calculate the residual energy generated by the time-dependent 
Hamiltonian. Finally we discuss possible experimental realizations of this work. 

PACS numbers: 64.60.Ht, 05.70.Jk, 64.70.Tg, 75.10.Jm 



I. INTRODUCTION 

Introduced in the 1930s independently by Landau and 
Zener [H, HJ, the Landau-Zener (LZ) transition formula 
provides an exact expression for the excitation probabil- 
ity in the final state when two levels approach each other 
due to a linear variation of the diagonal terms of a two- 
level Hamiltonian with an avoided level crossing at the 
point of minimum gap. Even after seventy years, the LZ 
formula is being applied extensively in problems over a 
wide range of modern physics, including neutrino oscil- 
lations, atomic physics, quantum optics, and mesoscopic 
systems Q. Recently, it has found several applications 
in quantum computations Q and adiabatic quantum dy- 
namics of many-body systems [fj. In particular, the gen- 
eral Kibble-Zurek (KZ) scaling [6-9] of the residual en- 
ergy or the defect density produced in the final state of a 
many-body Hamiltonian following a slow passage across a 
quantum critical point (QCP) [1(| has been established 
for a number of exactly solvable low-dimensional non- 
random spin models using the LZ formula 
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these systems factorize into decoupled 2x2 matrices in 
momentum space, and the dynamics gets mapped to a 
set of decoupled LZ problems. We note that there are 
several other studies of quenching dynamics of different 
critical systems [IB - [l7j . Quenching through a multicrit- 
ical point [HI , repeated passage through a QCP [l9[ and 
periodic variation of a parameter [2(| , quenching dynam- 
ics of a disordered spin chain [2l[ , and the dynamics of an 
open system coupled to a heat bath [22| have also been 
studied. 

In this paper, we address the following question: how 
does the transition probability of a LZ problem get al- 
tered if we introduce a waiting time t w at the minimum 
gap? Here, starting from the ground state at t = — oo, 
the system is brought to the point of minimum gap (2A) 
at time t = by linearly varying the diagonal terms of 
the Hamiltonian at a rate l/(2r). At the minimum gap, 



the system is allowed to evolve without any external driv- 
ing for a time t w , after which the linear variation is again 
resumed up to t = oo. We study how the probability of 
finding the system in the excited state at t = oo, given 
exactly by exp(— 2itA 2 t) for the conventional LZ prob- 
lem [1J , gets modified due to the additional time scale t w . 
We then extend our results for the Landau-Zener prob- 
lem with waiting to a many-body system whose quench 
dynamics through a QCP can be viewed as a set of decou- 
pled LZ problems; the system is brought from its initial 
ground state at t = — oo to the QCP at t = 0, where 
the Hamiltonian is not changed for a time t w , before it 
is again varied in time up to t = oo. As the relaxation 
time of the system diverges at the critical point, the re- 
sponse of the system to any perturbation becomes in- 
finitely slow; hence the system is no longer able to follow 
the ground state and therefore excitations are produced. 
As mentioned already, there have been several studies in 
this field recently, but the effect of waiting for a time t w 
at the critical point of a many-body system has not been 
studied. In a many-body system, it is more meaningful 
to look at quantities like the residual energy e r defined as 
the difference in energy between the actual state reached 
and the true ground state at the final time, or the density 
of defects (wrongly oriented spins) n which is obtained by 
integrating the probability of excitations over all modes. 

We will study the effect of waiting on the residual 
energy e r in the final state of the one-dimensional Ki- 
taev model following a quench through the QCP along 
with waiting. We study the possible correction to the 
KZ power law scaling n ~ i/ T du /^ z+1 \ where d is the 
spatial dimension and v and z are the correlation length 
and dynamical exponents associated with the QCP [1(| 
across which the system is swept. To the best of our 
knowledge, these questions have not been addressed be- 
fore from a theoretical point of view, although experi- 
mental results with waiting for single molecular magnets 
Mni2Ac are available [23|. The possibility of experimen- 
tal realizations of quenching dynamics with waiting in 
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optical lattices serves as another motivation. It is inter- 
esting to note that a similar concept is used in Ramsey 
spectroscopy where a molecule is subjected to an oscillat- 
ing perturbation for a time T\ which induces transitions 
between two specific levels of the molecule. The pertur- 
bation is then switched off for a time T2 after which it 
is again switched on for a time T\ and the probability of 
excitations is found (24|. Wc note that LZ sweeps have 
been used to generate coherent superpositions in quan- 
tum optical experiments [1^] . 

Our results can be summarized as follows. The effect 
of a waiting time t w can be understood by visualizing 
the dynamics in two parts: from t = —00 to and from 
t = to 00. If the solution of the first part is known, 
the second part can be solved by applying the idea of 
time reversal on the first part; in this way, given the 
amplitudes of the two basis states at t = 0, one can 
find the values of the same at t = 00. This reproduces 
the exact result for the conventional LZ problem. We 
then apply the method to a problem with waiting at the 
minimum gap to obtain an exact expression for the ex- 
citation probability in the final state; this gives simple 
forms for both the diabatic (A 2 r — > 0) and the adiabatic 
(A 2 r — > 00) limit. The probability of excitations exhibits 
a sinusoidal behavior in both cases but with different pre- 
factors and phases. The method is then used for a many- 
body system, namely, the one-dimensional Kitacv model 
[26j to obtain the residual energy. We quench this sys- 
tem through a QCP by linearly varying the anisotropy 
in the interaction, d J_ /dt = 1/r [13l . |27| , with a wait- 
ing time t w at the QCP. We show that for t w jy/wT <C I , 
the residual energy shows an exponential decay with t w 
given by e r ~ (a/^/r)[l + 6cxp(— ct w /y/r)]. We find that 
the parameters a, b and c obtained by an approximate 
analytical calculation give a good fit with the results ob- 
tained by numerically solving the Schrodinger equation. 

The outline of this paper is as follows. We describe our 
method of solving the Landau-Zener method in Sec. Ill Al 
and obtain an expression for the excitation probability in 
the presence of waiting in Sec. Ill Bl In Sec. IIIII we study 
the waiting problem for a many-body system taking the 
example of the onc-dimcnsional Kitaev model. 



II. LANDAU-ZENER PROBLEM 
A. Landau-Zener revisited 

To illustrate our method, let us revisit the conventional 
two-level LZ problem. Although the model has been ex- 
actly solved and the evolution matrix is known exactly 
for all times [2^, [2^] , and it has also been studied within a 
rotating wave approximation (30| . we describe a method 
below which can be easily generalized to study LZ dy- 
namics with waiting. 

The state ip{t) = C x (i)|l) + C 2 (t)|2), where |1) and |2) 
are the basis states (the initial and final ground states, 
respectively), evolves according to the Schrodinger equa- 
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where A is chosen to be real without any loss of general- 
ity. With the initial condition \C\{— oo)| 2 = 1, the wave 
function at t = is given by ip(t = 0) = a\l) + /3|2), 
where 
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with y = A 2 t/2 J2S, [3lJ . Henceforth, (•, -) T will denote 
the transpose of a given row vector. Since tf>(— 00) = 
(1,0) T evolves to i/j(0) = (a,f3) T , orthogonality implies 
that ^(-00) = (0, Vf must evolve to -0(0) = (/?*, -a*) T , 
up to a phase. Using properties of the Gamma functions 
|32l ]. one can show that |a| 2 + |/3| 2 = I as desired, and 
\a\ 2 -\P\ 2 = e-* A2r . 

Let us now ask: what wave functions at t = will 
evolve to (I,0) T and (0,I) T at t = 00? To answer this 
question, let us multiply the Hamiltonian in Eq. (JfJ by 
a z on both sides, which gives 



-t/(2r) A 
A t/(2r) 



(3) 



Clearly, by substituting t' = —t and ip'(t') = a z ?p(t) in 
Eq. (|3|), we recover Eq. (p}. Thus, the dynamics occur- 
ring in Eq. (Q]) from t — to —00 is the same as in Eq. §3§ 
with t' going from to 00 and tp replaced by a z tjj. Since 
ip(—oo) = (1, 0) T evolves to ip(0) = (a, (3) T , the above ar- 
gument shows that ip(0) = cr z (a, f3) T — (a, ~(3) T evolves 
to ip(oo) = (I,0) T through the Hamiltonian given in 
Eq. ([1}. By similar arguments, or by orthogonality, we 
see that ip(0) = a z ((3*,-a*) T = (f3*,a*) T evolves to 
^(00) = (0, 1) T , again up to a phase. We can now find 
the probability of ending in the excited state |1) at t = 00 
as follows. We can write ^(0) = (a,/3) T as 
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(MM/31 2 ) 



{a*, -13*) + 



/3* 



(4) 



where, in the first line, we introduced an identity opera- 
tor using the orthonormal basis (a, — (3) T and (/?*, a*) T . 
Since we now know the evolution of (a, —(3) T and 
(/?*, a*) T from t = to 00, we see that 



i<{oo) = (\af 
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(5) 



where there may be an unimportant phase difference be- 
tween the two terms. The probability of excitations at 
the final time is the probability to be in state |f) and is 
thus given by p = (|a| 2 — |/3| 2 ) = e - 27rA " r which is the 
exact expression for the LZ transition probability 
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B. Landau-Zener with waiting 

We now apply our method to the LZ problem with 
waiting. Within the time interval [0, t w ], the eigenvectors 
of H are given by (1, ±l) T /v / 2 with eigenvalues ±A. So 
the wave function changes from ip(0) = (a,/3) T to 



ip(t w ) 



a 



-»At„ 



a ~ P AAt w 

2 



1 

-1 



Inserting the identity operator used in Eq. (j4|), we get 



0* 



CM) 



ip{t v 



The state at t — oo can again be obtained by using the in- 
formation about the evolution of (a, —f3) T and (/?*, a*) T : 
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Hence the excitation probability is given by 

p tw = \{a*,~(3*)yj{t w )\ 2 

= [{\a\ 2 - |/? 2 |) cos(At U) ) - i(a*p - ap*)sm{M w )]' 



(6) 



This expression simplifies in the limits A 2 r — > and oo. 
For A 2 r = 0, we have pt w = cos 2 (At u) ); this is expected 
as the system does not get any time to evolve and there- 
fore remains in the state |1) up to t = 0, then oscillates 
between the states |1) and |2) from t = to t w , and 
then remains in the superposed state reached at t w for 
all t > t w . If A 2 r < 1, one gets i(a*/3 - (3* a) ~ Ix^/ttt 
which leads to an approximate expression 
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cos 2 [A(i„ 



7TT)]. 



(7) 



Comparison with the case A 2 r = shows a decrease in 
amplitude from 1 to e _7rA T and a phase shift of \Jttt. 



In the adiabatic limit A 2 r 



oo, a 



0, and 



one can use the asymptotic expansions of the Gamma 
function [32j to obtain the excitation probability 
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sin 2 (A< u) ). 



(8) 



In the next section, we consider waiting in a many-body 
problem, namely, the one-dimensional Kitacv model. 



III. WAITING IN KITAEV MODEL 

We now use the above results to study the quench- 
ing dynamics with waiting at the QCP of the one- 
dimensional Kitaev model given by the Hamiltonian 

MM 



(9) 



where j refers to the site index and a x,v are the Pauli 
matrices. The model can be exactly solved in terms of 
Jordan- Wigner fermions [33| defined as 



2j-l 2j 
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'2j+l- 



(10) 



Going to momentum space with %pk = (a,k,bk) T , where 
afc and bk are the Fourier transform of aj and bj, and 
performing an appropriate unitary transformation, the 
Hamiltonian decouples into a 2 x 2 form given by [27| 



H = 



tt/2 

E 

fc=0 



where Hk = 2 



H k ip k , 

J- sin(fc) J + cos(fc) 
J+ cos(fc) — J_ sin(fc) 



(11) 



J± = J\ ± Ji and k ranges from to 7r/2. The vanishing 
of the gap (= A^J J 2 sin 2 k + J 2 , cos 2 k) for the mode k = 
tt/2 at J_ = si gna ls a quantum phase transition of 



topological nature [34[, with v = z = 1. 

Setting J + = 1, we now apply the quenching scheme 



J-(t) = t/r for - oo < f < 
= for < t < t w 
= (t — t w )/r for t w <t< oo, 



(12) 



which incorporates waiting for a time t w at the QCP. The 
Schrodinger equation for each mode is 



J- (t) sin k 
cos k 



cos k 
-J- (t) sin k 



Mt), (13) 



so that the excitation probability for each mode is given 
by the modified LZ formula in Eq. (J6)). In Fig. [TJ we 
compare the exact analytical expression given in Eq. ((6]) 
with its corresponding approximate form for small A 2 t 
given in Eq. Q; the picture in the figure is reminiscent 
of Ramsey fringes. Noting that A 2 r = tcos 2 kj sinfc ~ 
0.00012 for the mode shown in Fig. [TJ we can use the 
diabatic limit result in Eq. (0, 



_ — 7TT COS 



fc/sinfc cos 2[2cOSfc(^ + - 



sin k 



)]■ (14) 



This gives a good understanding of the peak heights and 
phase shift in Fig. [T] 

The variation of Pk,t m vs k for different t w (shown in 
Fig. [2] with t = 10 ) shows secondary maxima, the peak 
heights of which increase as t w increases, in contrast to 
the conventional quenching case (27| . The increase in the 
number of maxima with increasing t w can be explained 
by the expression in Eq. (|14l) which has maxima at 



cos k 



2(t„ 



1FF/2Y 



m = 0,1,2, 



(15) 



in the limit of small k. With increasing t w , Eq. (|15|) is 
satisfied by more and more values of k which still satisfy 
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FIG. 1: Variation of probability of excitations pk with waiting 
time t w , for r = 1 and k = 1.56. The dots correspond to the 
exact solution given in Eq. I|6"]l. whereas the line corresponds 
to the approximate analytical expression given in Eq. (0. 
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FIG. 2: (Color online) Variation of probability of excitations 
Pfc,t„ with k for r = 10 obtained by numerically solving the 
Schrodinger equation, for t w = (solid red line), 4 (dotted 
black line) and 8 (dashed blue line). 



the condition that cos k is small enough so that the peak 
height given by the exponential pre-factor in Eq. (|14[) is 
not very small. Further, for a given value of m in Eq. 
(fT5|) . cosfc decreases as t w increases which justifies the 
increase in the peak height given by Eq. (|14p . 

The residual energy per site, e r , is given by the differ- 
ence of the expectation value of the operator 



O = 



'2m-\ u 2m 



(16) 



between the many-body state that is actually reached 
and the true ground state of H at t = oo (this is also the 
ground state of O); N denotes the number of sites. We 
find that e r is given by 



^ 2 dk 
2tt 



8 sin k pk. 



(17) 



Although this expression cannot be evaluated analyti- 
cally in general, one can obtain an approximate cxpres- 



FIG. 3: Variation of residual energy e r with waiting time t w , 
for r = 10. The line corresponds to the integral given in Eq. 
(fl7|) . while the dotted line corresponds to Eq. (fT8| . 



sion when r ^> 1 and t w / ^Jttt <C 1. Most of the contribu- 
tion to the integral in Eq. (fTTj) then comes from k close 
to 7r/2 in Eq. (fT4"]) . Approximating cosfc ~ ir/2 — k and 
sin k ~ 1 , redefining k = n/2 — k and finally extending 
the limits of integration from [0, 7r/2] to [0, oo], we obtain 



dk 



1 + cos[4fc(t w + -yfnr)] 



1 



0.32 



1 + g-^t^+iv^) 2 /^) 



[1 



0.37 e 



-2.3t™/v/T 



(18) 



where we have used the approximation t w /yjlrr 1 in 
the third line. In Fig. [3l the numerical results obtained 
by solving the Schrodinger equation arc compared with 
the approximate analytical expression given in Eq. (|18|) . 
This comparison shows that the numerical and analytical 
expressions are in good agreement. 

The decrease in the residual energy as a consequence of 
waiting can be understood as follows. For slow driving, 
only modes close to the critical modes contribute. For a 
mode with momentum fc, the frequency of oscillations of 
|Ci(£)| 2 between the two levels during waiting is propor- 
tional to Aft = 2cos/c; this vanishes as k approaches the 
critical value 7r/2, leading to a diverging time period T. 
Further, the vanishing of the off-diagonal element in the 
Hamiltonian in Eq. (flU)) implies that the modes close to 
fc = 7r/2 will remain close to the excited state under time 
evolution, i.e., \C\\ will remain close to 1. Now consider 
the variation of the diabatic excitation probability \C\\ 2 
between t = and t = t w . The waiting gives \C\\ the 
time to oscillate to | C 2 1 . Since we are considering only 
values of t w much smaller than the time period T, we en- 
counter only the decreasing part of the oscillating |Ci| 2 
such that when the variation of the parameter is again 
started at t w , \Ci{t w )\ 2 < \Ci(t = 0)| 2 . This reduces the 
probability of excitations at the final time. 

The advantage of studying the waiting problem in the 
Kitaev model is that the minimum gap for all the modes 
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occurs at the same time (t = 0) which makes the analyt- 
ical calculations easier. In many other models, the mini- 
mum gaps for the different modes, given by the vanishing 
of the diagonal term of the 2x2 Hamiltonians, do not 
occur at the same time. This makes it difficult to specify 
the precise time at which the waiting should be initiated 
so that our analytical results can be applied. But one 
can prove that in the limit of large t, the minima for all 
the modes approach the time at which the minima of the 
critical mode occurs. For large r, we therefore expect our 
analysis to go through with small corrections to the prob- 
ability of excitations given in Eq. (flU) . To understand 
why this is so, consider the quenching of the transverse 
magnetic field h in the transverse field XY model fl3 |. 
where the diagonal term of the equivalent 2x2 matrix 
is /;, + Jcosfc and the off-diagonal term is 7sinfc, with 
J = J x + J y and 7 = J x — J y . By expanding the di- 
agonal term about the critical mode k = 0, multiplying 
the corresponding Schrodinger equation with ypr and re- 
defining t 1 = t/ y/r, the diagonal term can be rewritten 
as t' + 1/(27 y^r), where the characteristic momentum 
scale is given by 1/ (7 s/r) by the LZ tunneling formula. 
Clearly, in the limit r — > 00, the minima for the mode 
occurs approximately at t' = 0. These arguments are 
applicable to the modes close to the critical mode; the 
modes far away from the critical mode anyway do not 
contribute to the residual energy or defect density as the 
excitation gap is very large. 

The waiting problem is also interesting from an exper- 
imental point of view. In Ref. l23l . the relaxation dynam- 
ics with waiting at a resonance was studied for single 
molecular magnets (SMM) called Mn^Ac, where each 
molecule has total spin S = 10. As the magnetic field in 
the ^-direction is varied from a large negative to a large 
positive value, S z changes from 10 to —10. In Ref. |23j |. 
the magnetic field is swept to a resonance value, start- 
ing from a field of —6 T, where it is held for different 
waiting times causing tunneling of the spins; eventually 
the field is brought back to its initial value. The number 
of molecules which have tunneled through the barrier, 
measured using electron paramagnetic resonance tech- 
niques, shows nearly a stretched exponential decay with 
t w , where the decay constant is governed by the relax- 
ation time of the system. Though there are limitations 
in mapping SMM to a LZ problem for all quenching rates 
[35[ | , it is worth noting that there have been experimental 
studies on the effect of waiting at the resonance, and im- 
portant quantities like the relaxation time of the system 
can be obtained from the waiting problem. Although, in 
our case, J_ is driven to a final value of +00 at t = +00, 
we do observe a qualitatively similar behavior, i.e., an 
increase in the tunneling probability to the second state 



with increasing waiting time. We believe that similar ex- 
periments with waiting and forward driving of the mag- 
netic field can be performed and it would be interesting 
to compare the results with our predictions. However, 
if the non-linear term of a SMM Hamiltonian [35| dom- 
inates, one may need to look at non-linear LZ problems 
[3?| with waiting. 

It may be mentioned here that quenching with waiting 
can be studied if the Kitacv model can be experimen- 
tally realized using cold atoms and molecules trapped in 
an optical lattice as proposed in Ref. l37l In this pro- 
posal, each of the couplings can be independently tuned 
using different microwave radiations. It is possible to 
investigate the evolution of the spatial correlation func- 
tion of the operator ib n a n+r , defined in Ref. [TH as a 
function of various parameters, where a n and b n denote 
some Majorana fcrmions operators. This spatial correla- 
tion function depends on pk which we have already ob- 
tained for the waiting case. Then the evolution of defect 
correlations can be detected by spatial noise correlation 
measurements as discussed in Ref. l38l. 



IV. CONCLUSION 

To summarize, the technique proposed here not only 
provides an exact result for the standard LZ problem 
but also enables us to estimate the excitation probabil- 
ity for the dynamics with waiting at the minimum gap 
and the residual energy for the related quenching dynam- 
ics of some many-body systems like the one-dimensional 
Kitaev model. We can derive simple expressions in some 
limiting situations and the approximate analytical results 
arc in excellent agreement with numerical results when 
tw/\/TTT « 1. The arguments leading up to Eq. (Tl8l) in- 
dicate that the KZ scaling law \lr du ^ uz+1 ^ will generally 
remain valid in the presence of waiting, except that the 
function multiplying the scaling term has a piece which 
decays with increasing t w . Finally, we have discussed 
some possibilities for experimentally testing our results. 
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